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Abstract 


An  analytical  study  is  perfonned  to  determine  die  dynamic  response,  natural  frequencies 
and  mode  sh^ies,  of  deep  C(»nposite  cylindrical  shells,  including  the  effect  of  through  the 
thickness  shear  straia  The  DSHELL  finite  element  program  is  used  to  predict  the  first 
four  natural  frequencies  and  the  results  are  compared  to  a  reference  using  the  Galerkin 
technique.  The  program  was  extended  to  problems  considering  simply  supported-free 
boundary  conditiems.  As  is  well  know,  the  free  boundary  is  rather  difficult  to  represent 
using  the  Galerkin  technique  tqiproach.  The  mode  sh^s  are  created  by  plotting  a 
surface-contour  plot  of  the  eigenvector  output  from  DSHELL. 

A  linear,  free  vibration  analysis  was  performed  on  two  ccxnposite  panels  in  which  damping 
effects  were  assumed  negligible.  The  analyzed  panels  are  made  of  graphite/epoxy  (Gr/Ep) 
material,  having  different  ply  orientations.  The  frrst  panel,  used  as  a  baseline  in  this  study, 
has  an  arclength,  longitudinal  length,  and  radius  of  curvature  of  12  inches,  11  inches,  and 
12  inches  respectively.  A  [0®/-45O/445®/90O]5  ply  orientation  under  the  siiiq)ly 
supported-free  boundary  condition  was  considered.  Comparisons  between  this  baseline 
panel,  using  DSHELL,  with  previous  holographic  experimentation  and  analytical  studies 
(  STAGSC-1  finite  element  program  was  previously  incorporated )  were  found  to 
correlate  well.  The  percent  difference  was  observed  to  be  less  tha  8%  in  all  cases. 

Also,  it  was  found  that  accurate  modeling  of  the  composite  shell  panel  was  dictated  by  the 
natural  frequencies  and  mode  shapes.  From  this  study,  four  nodes  have  to  be  positioned 
per  each  half  sine  wave  formed  in  the  eigenvectors. 


XI 


For  the  second  panel,  the  curvature  and  the  span  to  thickness  ratio  were  varied  in  order  to 
measure  effects  on  two  ply  orientations,  [0®/90®]5  and  [-45®/+450]s,  under  two  boundary 
conditions,  all  edges  siiiq)ly  supported  and  simply  supported-free.  The  result  showed  that, 
as  the  shell  becomes  deeper,  the  frequency  becomes  smaller. 

Also,  findings  show  that  as  curvature  increases,  the  natural  frequencies  for  both  laminates 
increases.  The  effect  of  increasing  the  thickness  was  more  evident  at  the  shallower  end  of 
the  composite  shell.  The  shear  effect  was  evident  for  small  values  of  span  u>  thickness 
ratios  and  large  curvature  under  the  two  boundary  conditions  studied.  The  percent 
difference  between  DSHELL  and  the  Galo'kin  technique  was  generally  less  than  10%  in 
all  cases  compared.  For  this  study,  DSHELL  proved  to  be  a  very  useful  engineering  tool. 
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An  analytical  study  is  peifonned  to  determine  the  dynamic  response,  natural  frequencies 
and  mode  shapes,  of  deep  composite  cylindrical  shells,  including  the  effect  of  through  the 
thickness  shear  strain.  The  DSHELL  finite  element  program  is  used  to  predict  the  first 
four  natural  frequencies  and  the  results  are  compared  to  a  reference  using  the  Galerkin 
technique.  The  program  was  extended  to  problems  considoing  simply  supported-free 
boundary  condidmis.  As  is  well  know,  die  free  boundary  is  radier  difficult  to  represent 
using  the  Galerkin  technique  i^iproach.  The  mode  shapes  are  created  by  plotting  a 
surface-contour  plot  of  the  eigenvector  output  fit»n  DSHELL. 

A  linear,  free  vibratitm  analysis  was  performed  (xi  two  composite  panels  in  which  damping 
effects  were  assumed  negligible.  The  analyzed  panels  are  made  of  graphite/epoxy  {OrlEp) 
material,  having  different  ply  orientations.  The  first  panel,  used  as  a  baseline  in  this  study, 
has  an  arclength,  longitudinal  length,  and  radius  cS  curvature  dl  12  inches,  1 1  inches,  and 
12  inches  respectively.  A  [0O/-45O/+45®/90®]s  ply  orientation  under  the  simply 
supported-free  boundary  condition  was  considoed.  Cmnparisons  between  this  baseline 
panel,  using  DSHELL,  with  previous  holographic  experimentation  and  analytical  studies 
( STAGSC-1  finite  element  program  was  previously  incorporated )  were  found  to 
cmrelate  well.  The  percent  difference  was  observed  to  be  less  tha  8%  in  all  cases. 

Also,  it  was  fcmikl  that  accurate  modeling  of  the  composite  shell  panel  was  dictated  by  the 
natural  frequencies  and  mode  shapes.  From  this  study,  four  nodes  have  to  be  positioned 
per  each  half  sine  wave  formed  in  the  eigenvectors. 
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Fcm:  the  second  panel,  the  curvature  and  the  span  to  thickness  ratio  were  varied  in  order  to 
measure  effects  on  two  ply  orientations,  [0®/90®]s  and  [-450/+45®]s,  under  two  boundary 
conditions,  all  edges  simply  supported  and  simply  supported-free.  The  result  showed  that, 
as  the  shell  becomes  deeper,  the  frequency  becomes  smaller. 

Also,  findings  show  that  as  curvature  increases,  the  natural  frequencies  for  both  laminates 
increases.  The  effect  of  increasing  the  thickness  was  more  evident  at  the  shallower  end  of 
the  composite  shell.  The  shear  effect  was  evident  for  small  values  of  span  to  thickness 
ratios  and  large  curvature  under  the  two  boundary  conditions  studied.  The  percent 
difference  between  DSHELL  and  the  Galerkin  technique  was  generally  less  than  10%  in 
all  cases  compared.  For  this  study,  DSHELL  proved  to  be  a  very  useful  engineering  tool. 
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EFFECTS  OF  THICKNESS  AND  CURVATURE  ON  THE  NATURAL 
FREQUENCIES  OF  CYLINDRICAL  COMPOSITE  SHELLS 


L  Introduction 


Ccmqxmte  materials  have  reached  great  acceptance  in  many  technological  fields.  Due  to 
their  superior  straigtfi-to<weight  and  stiffhess-to- weight  ratios  ccMopared  wid)  traditional 
alloys,  cooaposite  materials  are  ideal  for  aerospace  structure  applications. 

Composites  are  typically  formed  in  layers  with  the  high  sttengdi  properties  usually 
oriented  in  the  in-plane  directions.  Thus,  comporites  can  be  made  lighter  and  thinner  than 
(XHivoitional  isotrc^c  matoials.  Anotho' advantage  of  composites  is  that  a  laminate  can 
be  designed  to  satisfy  a  particular  case,  that  is,  the  layup  can  be  "fcHmed"  to  get  the 
strengdi  and  stiffness  required  in  a  specific  qrpUcation. 

Conqxwite  shdl  structures  have  been  widely  used  in  aeioq)ace  plications.  Wing  sidns 
and  fuselage  panels  are  some  examples  of  these  aj^lications.  They  are  also  called  thin- 
shell  structures.  The  excessive  vibration  in  these  structures  could  cause  fatigue  failure  or 
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shortea  tfadr  useful  life,  and  thus  dynamic  analysis  becomes  inqxntant  in  ord^  to  {xevent 
coU^rae  of  these  diin-sheli  structures.  Since  the  strength  varies  fitom  ply  to  ply  in  a 
huninate,  composite  materials  are  complicated  to  analyze,  therefore  numerical  techniques 
must  be  considered  in  order  to  solve  these  problems. 


1.1  Background. 

Much  woric  has  been  done  in  the  area  of  con^site  cylindrical  shells  over  the  last  few 
years.  It  has  long  been  known  diat  Classical  Plate  Theory,  based  on  the  Kirchhc^-Love 
hypothesis,  tends  to  produce  large  errors  when  dealing  with  composite  materials.  The 
Kirchhcrff-Love  hypothesis  assumes  that  straight  lines  normal  to  die  undeformed  shell 
midsurfBce  remain  straight  and  normal.  In  other  words,  transverse  shear  strains  are 
neglected,  resulting  in  overestimates  of  the  natural  finquendes. 

Reissner  was  the  first  b>  recognize  the  need  to  include  the  effects  of  transverse  shear 
effects;  Mindlin  followed,  and  added  the  effects  of  rotatory  inertia.  The  so  called 
Reissner-Mindlin  theory  assumes  that  while  cross  sections  remain  plane,  they  are  allowed 
to  rotate.  This  theory,  however,  does  not  satisfy  the  boundary  conditions  of  zero 
transverse  shear  on  die  top  and  bottrxn  surfaces  of  the  shell.  That  requires  application  of 
a  correction  facttjr,  and  is  comnKmIy  accepted. 

Reddy  [13,14]  assumed  that  the  displacements  of  the  midsurfaces  are  cubic  functions  of  z. 
This  leads  to  a  parabolic  distribution  fn*  the  transverse  shear  strain,  and  does  not  requite 
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a  correction  factor.  Linneman  and  Palazotto  [12]  used  this  ^proach  in  developing 
sdutions  f(Mr  the  natural  frequencies  for  synunetric  cylindrical  coiiq)osite  shells. 

The  dynamic  response  of  cylindrical  composite  panels  has  been  investigated,  here  at  AFTT, 
in  previous  studies.  Walley  [24]  studied  the  dynamic  response  of  a  curved 
Graphite/Epoxy  (Gr/Ep)  panel  clamped  all  around  with  centered  cutouts  ( 2"  x  2",  2"  x  4", 
4"  X  4" ).  Cyr  [4],  using  the  same  panels,  studied  the  effect  of  cutout  orientation  on  the 
natural  frequencies  and  mode  shapes  of  those  panels.  Levraea  [7],  using  a  finite  element 
code  called  STAGSC-1,  investigated  the  eccentricity  effects  of  square  cutouts  on  the 
natural  frequencies  and  mode  shtqres  of  curved  oxnposite  panels.  And  Silva  [16] 
investigated  a  con^x)site  cylindrical  shell  panel  under  transverse  load  with  through  the 
thickness  and  snapping.  Tsai  aiKl  Palazotto  [22]  investigated  the  non-linear  vibration  of 
cylindrical  shells  with  high-order  shear  deformation  dieory. 


1.2  Objectives. 

Using  certain  panel  geometries,  this  study  will  investigate  the  thickness  effects  cm  the 
natural  foequencies  and  mode  shtqies  for  different  boundary  conditions  and  different  ply 
(mentation. 

There  were  two  different  size  panels  ccmsidered  in  this  study,  the  Erst  erne  used  for 
con^arisem  had  an  arclength,  height  and  radius  of  curvature  of  0.305  m  (12  inches), 
0.2794  m  (1 1  inches)  and  0.305  m  (12  inches)  respectively.  One  ply  lay  up  was  analyzed, 
[0O/-45®/+45®/90O]s  in  which  each  ply  had  a  thickness  of  0.(X)0127  m  (0.005  inches). 
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giviiig  a  total  of  0.002032  m  (0.08  inches)  in  parcl  thickness.  The  second  panel  had  an 
arclength,  height  and  radius  of  curvature  of  0.508  m  (20  inches);  two  ply  orientations 
were  analyzed,  [0®/90®]s  and  [-45®/+45®]s.  Each  ply  was  0.(X)635  m  (0.25  inches) 
thick,  giving  a  total  panel  thickness  of  0.0254  m  (1  inch).In  order  to  facilitate  ply 
notation,  the  panels  mentioned  will  be  called  Panel  1,  Panel  2a  and  Panel  2b  respectively. 

The  boundary  conditions  considered  for  both  panels  were  Simply  Suppmted  along  the 
vertical  edges  and  free  along  the  horizontal  edges  (SFSF),  and  Simply  Supported  all 
Around  (SS).  The  DSHELL  finite  elemoit  code  was  used  to  predict  the  Erst  four  natural 
frequencies.  The  dynamic  capability  df  this  program,  in  predicting  the  natural  frequmcies 
and  mode  shapes,  was  evaluated  conqnuing  to  a  reference  using  the  Galeikin 
technique. 


L3  Analytical  Approach. 

The  numerical  analysis  was  conducted  using  the  DSHELL  ccnnputer  code.  Previous 
analyses  have  been  successfully  conducted  at  AFTT  using  this  finite  element  program  but 
never  a  dynamic  frequency  respmse  of  cylindrical  composite  shells. 

DSHELL  assumes  a  parabolic  shear  strain  which  varies  through  die  shell  thickness  and 
vanishes  at  die  tq)  and  botumi  surface.  Although  classical  plate  and  shell  theory  for  thin 
isotropic  structures  ignore  shear  stress  through  the  thickness,  it  is  not  appropriate  to  do  so 
with  composite  shells.  The  coupling  of  extensitnial,  bending,  and  shear  strain  must  be 
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taken  into  consideration. 


DSHELL  contains  a  built-in  grid  generator  which  was  used  to  generate  all  of  the  finite 
element  meshes  required  for  this  investigation.  This  program  was  selected  fm*  use  in  this 
stixly  to  continue  verification  of  its  dynamic  capabilities. 
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n.  Theory 


hfost  aenrapace  structures  are  considered  continuous  systems  and  have  distributed 
material  properties.  An  aircraft  wing  can  be  considered  as  a  beam  with  nonuniform  mass 
per  unit  length  and  flexural  stiffness  while  a  fuselage  panel  can  be  mealed  as  a  cylindrical 
shell.  Probkms  concerning  the  vibration  of  shells  are  considerably  more  complicated  than 
their  counterparts  for  beams  w  plates.  Primarily,  this  is  caused  by  the  effects  of  the 
curvature  on  the  shell  equations  and  thus  on  the  dynamic  behavior.  For  shells,  membrane 
and  flexural  deformations  are  coupled,  and  any  themy  must  consider  these  effects 
simultaneously. 


2.1  Natural  Freqiiendes  of  Shdb 

Li  the  classical  laminated  shell  themy,  through  the  thickness  shear  defoim^n  is 
neglected  according  to  the  Kirchhoff-Love  hypotesis  that  plane  cross  sections  remain 
plane  before  and  after  deformatim.  Mindlin  introduced  a  theory  that  allowed  the  cross 
section  remain  plane  but  rotate  from  the  normal  with  respect  to  the  mid-surface  afto* 
deformation.  The  assunqition  of  no  cross  sectional  warping  introduces  error,  especially  at 
the  top  and  bottom  surfaces,  since  the  boundary  condition  of  zero  transverse  shear  is  not 
matched. 

Reddy  [13]  introduced  a  new  theory  in  which  the  displacements  of  the  mid-surface  were 
cubic  fuiKtios  of  z.  This  lead  to  a  paraboUc  transverse  shear  distribution  wherein  the 
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stfains  are  maximum  at  the  shell  mid-surface  and  are  zero  at  the  tq)  and  bottom  surfaces, 
satisfying  the  required  thickness  boundary  condition.  Fig  2-1,  shows  the  transverse  shear 
concq>tsdisscussed  above. 

BEFORE  DEFORMATION 

Plane  Cross  Section 


AFTER  DEFORMATION 
Kiichhoff  Theory  _ 


Mindlin  Theory 


Parabolic  Themy 


Fig.  2.1  Transverse  Shear  Strain  Theories 


To  achieve  die  desired  parabolic  shear  distributicMi,  a  higher  ordo*  displacement  field  is 
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required.  This  theory  has  been  studied  by  Reddy  [13]  and  also  by  Palazotto  and  Dennis 
[11];  both  have  presented  higher  order  shear  theories  incorporating  transverse  shear 
effects,  and  this  is  certainly  the  most  inqwrtant  feature  of  DSHELL.  Tsai  and  Palazotto 
[22]  incorporated  the  capability  to  performed  dynamic  analysis  into  this  approach. 

A  description  of  the  main  features  involved  in  DSHELL  will  be  explained  in  order  to 
clarify  the  theory  behind  this  finite  element  program. 


2.2  The  K  and  M  matrixes 

This  section  is  directed  to  the  development  of  the  appropriate  vibration  equations. 

The  governing  equation  for  the  dynamic  behavior  of  a  continuous  system,  derived  from 
the  Hamilton's  principle  is  [10:44] 

m  Xit)  +  c  Xit)  +  k  X{t)  =  F{t)  (2.1 ) 


where 

m  =  system  mass 

X(t)  =  displacement 

c  =  system  damping 

X(t)  =  velocity 

k  =  system  stiffness 

X(t)  =  acceleration 

t  =  time 

F(t)  =  external  forcing  function 

Fcm*  free  vibration  there  are  no  applied  forces  and  therefore  F  (t)  =  0.  Structural  damping 
is  based  on  the  energy  losses  caused  by  internal  friction  as  the  panel  deforms  and  is 
proportional  to  the  amplitude  of  displacement  [10:402].  The  extremelly  small 
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displacements  that  take  place  in  these  cases,  compared  with  the  large  stiffness  of  the 
system  [21],  justified  neglect’ig  the  damping  effects,  therefore  c  =  0.  Thus,  equation 
( 2.1 )  becomes 


m  X(t)  +  k  X(t)  =  0  continuous  system  ( 2.2 ) 


The  [K]  matrix,  for  an  element  modeling  a  composite  material  panel  considering  large 
displacements,  is  given  by  [22] 


K  (u)  =  Xo  + 


ic,(u)  ^  K,(ur 
2  3 


(2.3) 


where  Ko  is  a  constant  stiffness  matrix,  ilTi  is  a  stiffness  matrix  related  to  linear 
displacement,  and  X2  is  a  stiffness  matrix  related  to  quadratic  displacement.  For  this 

analysis,  the  Xo  matrix  will  be  the  only  component  used.  The  consistent  mass  matrix  is 
obtained  from  [22] 


^  =  I  ll!  P  (2.4) 

whereCft-i  and  are  positions  at  the  top  and  bo^m  of  the  ^  th  layer,  £2  is  the  domain  of 
the  mid  surface  and  L  denotes  the  number  of  layers  in  the  laminate.  [/?]  is  a  matrix  that  is 
a  function  of  ^  only.  A  full  discussion  this  matrix  is  referred  to  reference  [22],  where 
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[A^]  is  the  shape  function  matrix  and  p  is  mass  per  unit  volume  (density).  The  consistent 
mass  formulation  results  in  a  fully  populated,  positive  defmite  mass  matrix. 

The  finite  ehnent  method  requires  the  discretization  of  a  continuous  system  into  small 
elements  in  which  the  deformed  geometry  is  represented  by  displacements  at  selected 
locations  (nodes)  within  the  body.  Equation  ( 2.2 )  can  be  written  as 

[Af]  (d)  +  [AT]  [d]  =  0  discrete  system  (2.5) 


where 

[M]  =  mass  matrix 

[d]  =  nodal  displacement  vector 

[K]  =  stiffness  matrix 

•  •  »• 

(d)  =  nodal  acceleration  vector 

Fot  harmcxiic  vibration  it  is  assumed  that  all  nodal  DOF  are  in  phase  and  that  ( d  }in 
Equation  2.5  is  given  by  [3:308] 


[d]  =  (2.6) 

where  d^  is  the  displacement  amplitude,  <o  is  circular  natural  frequency,  and  t  is  time. 
Substituting  Equation  2.6  into  Equation  2.5  gives 

[K]{d^}  [do)  =  0  (2.7) 

Rearranging  and  changing  notation,  one  obtains 
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[K-X,M]{^,]  =  0 


(2.8) 


where  Xi  =  od,-^  ((rad/sec)^)  and  the  brackets  have  been  dropped  off  of  K  and  M  for 
convenience  (AT  and  Af  are  still  matrices).  Fm*  linear  vibration  analysis,  AT  and  Af  are 
independoit  of  displacement  and  frequency.  X^  is  the  eigenvalue  and  (ijkjisthe 
eigenvector  (mode  shape)  for  the  i^  mode.  Each  eigenvector  will  represent  a  unique 
shiqw  but  with  arbitrary  amplitude. 

Equation  2.8  rq»esents  the  generalized  eigenvalue  problem  for  linear  fifce  undamped 
vibration  [1:517].  If  AT  and  Af  are  (n  x  n)  matrices,  then  the  solution  to  Equation  2.8 

yieldsXj  and  for  i  =  1, 2, 3,  ...n.  The  eigenvectors  are  independent  and  thus  form  a 
basis  for  the  eigenspace  [1 :49].  The  eigenvectors  can  be  airan^  into  a(nxn)  modal 
matrix  given  by 


^  =  [  («,}  {<l>2)  ...{9J]  (2.9) 

where  the  columns  of  the  modal  matrix  are  the  eigenvectors  of  the  system.  This  allows 
Equatitm  2.8  to  be  written  as  [1:558] 

ATO)  =  AfOA  (2.10) 

Here,  A  is  a  diagonal  matrix  of  the  eigenvalues,  Xg  and  their  location  along  the  diagonal  is 
determined  by  Xj,  <  X22  <  ...  <  X„  .  The  location  of  { <!>, }  within  O  must  be  consistent 
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with  die  locaticMi  of  in  A;  that  is,  { <|>j }  is  the  conresptHiding  vector  for  X,, . 


To  solve  the  eigenvalue  problems  DSHELL  uses  a  subspace  iteraticm  scheme.  The 
procedure  iiiq)lemented  in  DSHELL  is  based  on  the  subspace  iteration  method  developed 
by  Bathe  [1:672-682]. 


23  Subspace  Iteration  Method. 

The  inirpose  of  this  section  is  tt>  discuss  the  subspace  method  incorporated  in  the 
DSHELL  code  to  find  the  eigenvalue-ei^nvector.  Appendix  Aindicates  some  additional 
details. 

The  basic  objective  in  the  subspace  iteration  method  is  to  solve  for  the  lowest  p 
eigenvalues  and  corresponding  eigenvectors  satisfying  [1:672] 

K<b  =  M<b\  (2.11) 

where  A  =  diag  (A,)  and  O  =  [<|>i,<|>2....,^,,  ] 

In  addition  to  Equation  ( 2.1 1 ),  the  eigenvectors  also  satisfy  the  orthogonality  conditions 

=  A;  =  /  (2.12) 

t^iere  /  is  an  unit  matrix  of  order  q  because  <b  stores  only  q  eigenvectors.  The  first 
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rdatioa  is  a  necessary  and  sufficient  conditkw  f(x  die  vectors  in  O  to  be  eigenvectors,  but 
die  eigenvector  orthogonality  conditions  in  the  sectHid  expression  are  necessary  but 
not  sufficient 

The  number  (tf  modes  defied  are  "p" ,  by  using  sub^ace  iteration  die  reduced  number 
of  "p"  eigenvectors  are  simultaneously  (ditained.  The  procedure  modifies  Bathe's 
subqiace  method  in  (Hdor  to  enhance  convergence. 

Subspace  iteration  reduces  an  ( n  x  n )  eigensystem  into  a  ( n  x  9)  system  where  "q"  is  as 
defined  before.  The  subspace  refers  to  the  space  spanned  by  the  reduced  set  of 
eigenvectors  space  [  1 :67S].  If  the  subspace  which  is  spanned  1^  the  reduced  eigenvectors 

is  denoted  ,  then  the  starting  iteration  vector  can  be  thought  of  as  spanning  the  £1 

subspace.  As  iterations  are  petfornKid,  the  E^  ( ksl,  2, 3,...,  q)  space  approaches  the  E^ 

space.  The  £,  qiace  is  in  effect  spanned  when  the  desired  accuracy  has  been  obtained. 

The  following  gives  an  outline  of  the  eigensoludon  technique  used  by  DSHELL.  Hrst,  a 
preliminary  set  of  independent  vectors  is  chosen  which  spans  a  subspace  of  the  complete 
set  of  dgenvectors  determined  by  the  actual  number  of  vectors  desired  ( i.e.  2p  or  p  +  8 
which  ever  is  smaller  ).  Next,  an  inverse  power  method  ( inverse  iteraticm )  with  shifting  is 
performed  ( the  shifting  is  used  to  overcome  any  matrix  singularities ).  This  transforms 

the  El,  space  into  space.  Now  the  projections  of  K  and  M  are  found  on  the  new 

reduced  eigenspace  ( subspace  ^*+1  )  yielding  K^^.,  and  respectively  [7].  The 
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reduced  (nxq)  ei^nvalue  problon  is  solved  using  Householder's  Method  with  a  IK 
algcxidun.  Finally,  a  check  for  convergence  of  the  eigenvalues  is  made  and  the  process 
iq)eated  until  adequate  ctmveigence  has  been  obtained.  The  criteria  for  determing 
ccmvergence  is  based  upon  the  ”k  +  1"  iteration  eigenvalues  being  ccmqiared  to  die  "k  " 
iteration  eigenvalues.  If  the  relative  error  is  below  a  predetermined  diieshold  then 
convei^nce  is  assumed. 


2A  Thickiiess  and  Curvature  Parameters. 

In  considering  a  cylindrical  shell,  thoe  are  two  geometric  parameters  by  which  a  shell  can 
be  described.  The  ratios  S/c  and  8/h  ate  used  to  explain  the  deepness  and  thickness  of  the 
shell  respectively  in  subsequoit  sections.  The  fcdlowing  figure  shows  a  typical  shell 
element  in  which  these  parameters  are  given. 


Fig  2.2  Typical  Shell  Element 

From  fig.  2.2,  5  is  the  (rffset  distance  measured  fircnn  the  cord  line  to  the  datum  surface  of 
die  shell,  h  is  the  diickness  of  the  shell  and  c  is  half  the  distance  of  the  cord. 
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Theconoeptofdeeimessandshaliownessisvery  useful  when  dealing  with  shells.  A  shell 
is  conskkied  shallow  when  the  ratio  S/c  is  small,  usually  S/c  ^  1/4,  [1 1;2Q2].  A  shell  with 
Wc  ^  1/4  is  considered  deep.  A  deep  shell  has  more  coupling  between  displacements 
functkms.  The  in-plane  displacements  u  and  v  are  affecting  the  rotation  about  a  normal 
to  the  shell  surface.  DSHELL  includes  these  effects,  and  thus  it  is  possible  to  determine 
wha[ithey  are  of  importance  by  studying  eigoivector  relations.  Inacklition,  a  shell  may 
be  considered  thick  when  the  ratio  S/h  <  1  [11:203]. 

According  to  the  objectives  of  this  study,  the  effects  of  thickness  and  curvature  tm 
composite  cylindrical  panels  will  be  evaluated,  therefore  the  following  figure,  considering 
the  case  of  the  [0^/90°]s  ply  orientation,  with  b/h  »  20  under  the  simply  supported  all 
around  boundary  condition  (discussed  in  Chapter  IV),  wiU  help  to  visualize  when  a  shell  is 
considered  thick  and  or  deep. 
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Hg  2-3  Typical  bdiavior  of  die  Fundamental  Frequoscy  for  the  [0^/90*^]$  ply 
orientation.  Simply  Supported  all  around  Boundary  CiMidition. 

The  ftdlowing  taUe  shows  the  values  for  d/c  luid  S/h  ratk>s  for  the  interval  b/h  and  the 
curvatures  showed  in  the  figure  above.  The  gemnetry  of  the  shell  and  the  parameters 
S/c  and  Mi  invcdved  are  given  in  figure  2-2. 


Rs20 

Osis/R 

5  =  R(l-cosfl/2) 

c 

5/c 

Ml 

b/h  =10 

28.640 

0.62 

4.95 

0.125 

0.621 

b/h»20 

57.290 

2.45 

9.58 

0.254 

2.44 

85.940 

5.36 

13.62 

0.393 

5.36 

b/h»40 

114.590 

9.19 

16.83 

0.546 

9.19 

b/h  =  50 

143.240 

13.70 

18.97 

0.722 

13.70 

1  n 

liBkl 

6  =  R(l-cos0/2) 

c 

5/c 

Ml 

b/h=10 

2.860 

0.062 

4.97 

0.0125 

0.062 

b/h  =  20 

5.730 

0.249 

9.97 

0.0252 

0.249 

b/h  =  30 

8.590 

0.562 

14.98 

0.0375 

0.562 

b/h  =  40 

11.460 

0.999 

19.96 

0.050 

0.999 

b/h  =  50 

14.320 

1.427 

23.84 

0.0598 

1.427 

Table  2-1.  Geometric  Parameters  of  the  Conqwsite  Shells 
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Accofding  to  this  table,  for  a  radius  R  =  20,  a  shell  with  a  ^Mun  equal  to  ten  (b/h  »  10)  is 
shallow  (8/c  »  0.125),  but  also,  at  the  same  time,  this  shell  can  be  con^dered  thick  since 
£^  =  0.62.  As  the  span  increases,  the  shell  becomes  deeper  since  S/c  is  greater  than  0.25 
but  also  it  becomes  diinner.  For  a  radius  of  curvature  R  =  200,  no  matter  what  the  span 
b,  the  shell  is  called  shallow.  It  can  be  considered  thin  cmly  for  values  b/h  ^  40. 

The  curvature  also  has  inq)ortance  in  the  behavior  of  die  frequency  response  for  the  case 
of  cylindrical  shell  panels.  According  to  fig.  2-3,  as  the  radius  of  curvature  decreases  to 
h/R  =  1/20,  the  frequency  increases.  This  behavio*  is  due  to  the  fact  that  the  in-plane 
displacements,  u  and  v,  have  more  effect  on  the  natural  finequencies  when  the  panel 
becomes  deeper  (when  the  radius  of  curvature  decreases).  Then  for  the  ratio  h/R  = 
1/200,  diese  effects  are  practicaUy  negligible.  In  c»der  to  undmtand  this  statement,  the 
next  figures  show  the  magnitude  of  the  in-plane  displacements  for  the  case  of  the 
[-45^/'M50]s  laminate,  sinqtly  supported-fiee  boundary  crnidition  and  b/h  =  20  for  two 
radius  of  curvature:  1/20  and  1/200. 


u  displacement  v  displacement  w  di^lacemoit 


Hg.  2-4  In-Plane  Displaconents  fm*  tiie  [0^/90^]$  laminate.  Simply 
Siq)pQrted-Free  Boundary  Condition,  b/h  =  20,  h/R  =  1/20. 
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udisplaconent  v  di^lacement  w  di^lacanmt 

Fig.  2-5  In-Plane  Di^lacemoits  fm*  the  [(yV90^]s  laminate,  Simply 
Supported-Ree  Boundaiy  Conditicm,  =  20,  h/R  =  1/200. 

From  the  ou^t  of  DSHELL,  it  was  found  that  for  h/R  =  1/20,  =  0.0522 

and  =  0.246  ,  while  for  the  ratio  h/R  =  1/200,  =  6.00522  and 

s  0.0245  .  It  means  that  die  in-plane  displacements  contribute  ten  times  more 

to  die  ovoall  eigenvector-eigenvalue  solution  for  the  h/R  -  1/20  than  the  ratio  h/R  = 
1/200.  It  should  be  noted  that  the  DSHELL  program  is  always  developing  eigenvectors 
and  eigenvalues  in  which  all  the  degrees  of  freedom  are  considoed,  as  repose  to  analyses 
in  which  in-plane  displacement  and  rotatmy  inertias  are  suppressed.  Thus,  by  looking  at 
die  maximun  mdinate  of  the  re^iectivedispiacements  functions,  one  can  obtain  an 
iqipreciation  of  the  ovoall  contribution  to  the  global  eigenvecttHr.  The  same  thing  was 
also  observed  fen*  the  rotation  contribution  and  Yj* 
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nL  Finite  Eleiiient  Analjfsis  and  Modelling 


3.1  Fanel  Properties,  Geometry  and  Aflsumptions. 

Hie  material  propefties  ami  geometry  of  the  panels  are  explained  in  this  section.  The 
analyzed  panels  are  made  of  graphite/epmy  (Gr/Ep)  material,  they  are  cylindrical  in  shape 
and  have  different  ply  orientation.  For  consistency  widi  previous  works  [7,8]aquasi- 
isotrcqiic  [0®/-45®/445®/90®]s ,  [0®/90®]s ,  ami  [-45®/+45®]s  ply  layup  were  selected. 
They  will  be  referred  as  Panel  1,  Panel  2a  and  Panel  2b  respectively.  Table  3-1  lists  the 
material  properties  for  the  panels  refnred  above;  the  properties  for  Panel  #  1  were 
obtained  fiom  experimentation  [7]  while  the  properties  for  Panel  #  2  were  taken  as  they 
were  assumed  using  the  Galeikin  technique. 


Mechanical  Properties 

Panel  #  1 

Panel  #2  (aandb) 

E,  [GPaj 

141.06 

144.79 

Ej  [GPa] 

9.238 

9.652 

[GPa] 

5.955 

4.136 

Gu  [GPa] 

5.955 

4.136 

Ga  [GPa] 

2.977 

2.068 

Va 

0.3131 

0.300 

V21 

0.0205 

0.020 

PlKg/m’l 

1612.02 

1522.3 

Table  3. 1  Material  properties  of  Composite  Panels 
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In  Panel  1  each  ply  iim  1.27  x  10*^  m  thick  resulting  in  a  panel  thickness  of  0.0Q2  m. 
Panels  2a  and  2b  had  a  ply  thickness  of  0.00635  m,  resulting  in  a  panel  thickness  of 
0.025401. 

The  curvilinear  orthogcMud  coordinate  syston  and  nooienclatuie  used  in  this  formulatitm 
of  die  laminated  cylindrical  shell  is  shown  in  Fig.  3-1.  The  x-axis  lies  along  the  straight 
dimension  of  the  panel;  the  s-axis  follows  the  circumference,  and  the  z-axis  is  everywhere 
nonnal  to  the  shell  middle  surface,  positive  toward  the  center  of  curvature.  The  surface 
formed  by  the  x  and  s  axes  lies  in  the  center  of  the  thickness  of  the  panel,  so  the  thickness 
comdinate  is  negative  on  the  outa  surface  and  positive  on  the  inner  surface. 
Displacements  along  the  x,  s  and  z  axes  are  u,  v  and  w  respectively.  The  angle  ^lecifies 
the  orientatitm  angle  of  eac!i  piy  in  the  laminate.  The  effective  panel  dimensions  are 
defined  as  the  dimensious  of  that  portion  of  the  panel  free  to  vibrate. 


Fig  3-1.  Shell  Panel  Dimensions  and  Coordinate  System 
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Two  boiuidary  comfitions  were  analyzed  for  both  panels:  Simple  Supported-Free  ( SFSF; 
circumferential  edges  are  free )  and  Simply  Supported  all  around  ( SS ).  The  panel 
chanusteristics  cited  in  this  section  were  used  to  create  part  of  the  finite  element  models 
( input  data  ). 


DSHELL  Program 

According  with  die  finite  element  analysis,  the  moddling  cylindrical  shells  have  been 
(kme  through  the  following  theories: 

a 

a )  Flat  Elonents:  Bane  bmding  elements  with  plane  membrane  elements 

b)  Curved  Elemoits:  Classical  shell  theory 

c )  Mindlin  type  elements:  special  forms  of  solid  elements 

The  DSHELL  finite  element  program,  used  in  the  present  study,  uses  the  type  of  elements 
described  in  part  (c).  This  finite  elemoit  code  was  originally  written  by  S.  Dennis  [1 1]  and 
modified  by  C.  T.  Tsai  and  Palazotto  [22],  The  uniqueness  of  the  DSHELL  program  is 
that  shear  strains  are  assumed  to  vary  parabolically  through  the  thickness  of  the  shell, 
vanishing  at  the  top  and  bottom  surface,  large  displacements  and  rotations  are  considered. 

DSHELL  also  incmporates  material  linearity.  In  structural  applications,  materials  are 
normally  restricted  to  the  linear  elastic  region  [6].  This  inq)lies  that  any  fiber  Ineakage 
during  large  displacements  and  rotations  is  ignored. 
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The  majority  ci  past  efforts  has  dealt  with  static  analysis  or  with  the  dynamic  analysis  but 
without  taking  account  of  the  eigenvectOT  problem.  Therefcne,  this  theas  will  serve  to 
further  validate  die  dynamic  modeling  capatnlides  of  the  DSHELL  program.  A  linear, 
free  vibration  analysis  was  conducted  on  the  Gr/Ep  circular  cylindrical  panels.  For  free 
vibration  analysis,  DSHELL  outputs  the  desired  natural  frequencies  and  ccnresponding 
eigenvectors.  They  represent  the  mode  shapes  of  the  panel. 

The  user  has  the  capability  of  selecting  the  type  of  element  under  study :  plate  or 
cylindrical  shell,  number  of  nodes  per  element  and  type  of  analysis;  ncmlinear,  linear  or 
eigenvalue  problem.  Also  the  number  of  natural  frequencies  required  and  type  of 
laminate  can  be  set  In  addition,  the  built-in  mesh  generatOT  was  used  to  establish  the 
finite  element  grid  for  all  the  panels. 

The  user  can  determine  the  number  of  elements  in  the  x  and  y  direction,  the  grid  spacing  in 
the  same  coordinate  directions  and  DSHELL  determines  the  nodal  coordinates  and 
(XMmectivity  matrix.  Any  of  the  degrees  d  freedom  at  each  node  can  be  restrained 
according  to  the  boundary  conditions  desired.  Finally,  the  material  properties  must  be 
input:  density.  Young's  modulis,  shear  modulus,  Poisson's  ratio,  number  of  plies,  ply 
orientation  angles,  ply  thickness  and  circular  cylinder  radius.  For  the  type  of  analysis 
perfrumed  and  the  panel  geometries  investigated,  DSHELL  proved  a  straight  forward 
implementation. 
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S3  EhHnent  Selection. 


Hie  DSHQX  finite  elrawnt  prc^ram  uses  a  36  DOF  element  developed  by  Dennis  [11] 
and  pictured  in  Rg  3-2.  The  seven  displacements  functions  m,v,  w,w,,  ,w,2  .'Pj  and  'Pj 
are  considraed  in  each  element.  The  mid-side  nodes  have  only  the  two  inplane  degrees  of 
freedom  u  and  v .  Lagrangian  bilinear  interpolation  is  used  for  4',  and  quadratic 
interpolation  is  used  fOTU  and  V,  while  nonconforming  Hermitian  interpolation  is  used 
for  w,  w,,  and  •  Although  a  shell  is  thr^  dimensional,  a  two  dimensional  formulation 
is  used  to  model  a  thin  shell. 


Fig  3-2.  36  DOF  Isoparametric  shell  element 
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where  :  ^  =  xfa  andt]  =  f /h,  the  natural  curvilinear  cotndinates  of  the  surface 


u  is  the  displacement  in  the  x  coordinate  direction 
V  is  the  displacement  in  the  s  coordinate  direction 
w  is  the  transverse  displacement  in  the  z  coordinate  direction 
w,,  is  the  slope  dwidx  at  the  node 
w,2  is  the  sl(^  dw/ds  at  the  node 

'P,  is  the  rotation  of  the  normal  to  the  shell  midsurface  in  the  x  direction 
( about  the  s  axis  ) 

is  the  rotation  of  the  normal  to  the  shell  midsurface  in  the  s  direction 
( abou'  the  x  axis ) 

The  overall  equation  defining  the  isoparametric  discretizaticm  of  the  continuum 
displacements  into  discrete  displacements  is  obtained  by  [19] 

|«|  =  [iV]  {q} 

(7x1)  ( 7  X  36 )(  36  X  1 ) 


3.4  Convergence  Study. 

A  convergence  study  was  carried  out  for  the  dynamic  response  of  composite  panels. 
This  study  was  conducted  using  the  36  DOF  element  described  before  and  Panel  #  1 


3-6 


(explained  in  section  3.1).  This  panel  had  been  analyzed  before  [7],  using  the  STAGSC-l 
finite  eleotent  code,  and  experimentally  by  Holographic  Interferomerty.  Since  results 
fitm  these  methods  were  available,  this  panel,  under  its  respective  boundary  condidm, 
was  chosen  <mly  to  prove  the  precisicm  of  the  DSHELL  code  in  predicting  the  natural 
frequencies  and  mode  shapes. 

The  DSHELL  built-in  grid  generator  was  used  to  create  all  element  meshes.  This 
resulted  in  meshes  composed  of  rectangular  elements.  The  following  pattern  was  used: 
the  number  of  elements  or  nodal  divisions  in  the  s  direction  (according  to  fig.  3-1)  was 
varied  while  the  number  of  elements  in  the  x  axis  kept  constant  Then,  when  the  number 
of  elements  in  the  5  directitm  were  kept  constant,  the  number  of  elements  in  the  x  axis 
were  varied. 

Bg  3-3  shows  the  behavior  of  the  firequoicies  wh«i  the  number  of  divisitms  in  the  x  axis 
were  10,  IS  and  30  elements  respectively.  It  is  seen  fix>m  these  curves  that,  as  the 
number  of  elements  increase  the  difference  between  reference  [7]  and  this  analysis 
becomes  smaller. 

Bg  3-4  show  the  behavior  of  the  firequencies  when  the  divisions  in  the  s  direction  were 
varied,  the  same  intervals  were  analynd  and  now  the  first  two  modes  were  affected  when 
more  than  15  divisions  were  used  in  that  direction.  The  higher  modes  ccmverged  to  the 
reference  soluticm. 


3-7 


DlttMWie* 

P«fo«nt 


Hg  3-3  Percent  Difference  in  the  Frequencies  by  Varying  the  x  Axis 


4S 

40 

35 

30 

DWwwie* 
P«rc«nt  20 

IS 

10 

6 

0 


10  15  ZO 

Nuntow  ol  Division*  in  tlw  s  axis 


Fig  3-4  Percent  Difference  in  the  Frequencies  by  Varying  the  s  axis 


According  to  this  analysis,  ten  divisions  should  be  used  in  the  x  axis  while  fifteen  elements 
were  required  in  the  s  direction.  When  these  results  woe  combined  in  a  mesh,  the 
solution  converged  to  the  reference  study  and  the  average  difference  was  about  4%. 
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Fig.  3-5  Typical  10  x  15  Element  Mesh 


Some  other  mesh  combinaticms  woe  used  and  they  are  included  in  the  next  table. 


Mode 

Number 

111111 

Mesh 

8x20 

Mesh 

10x15 

Mesh 
llx  12 

Mesh 

5  X  40 

1 

108.70 

95.36 

102.83 

104.8 

91.62 

2 

203.3 

208.42 

199.92 

219.3 

181.67 

3 

235.1 

248.21 

252.60 

268.2 

259.64 

4 

369.4 

368.37 

370.24 

398.2 

389.21 

Table  3-2  Natural  Frequencies  (Hz)  from  Convergence  Study  for  Panel  #  1 


Hgures  3-6  and  3-7  show  the  ixiode  shapes  obtained  for  two  meshes  analyzed:  the  mesh 
used  and  the  zeference. 


Third  Mode  Shape 


Fourth  Mode  Shape 


Fig  3-6.  Mode  Shapes  for  Panel  #  1  according  to  the  reference 
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First  Mode  Shiqie 


:•  Second  Mode  Shape 


Third  Mode  Shq)e 


Fourth  Mode  Shq)e 


Fig  3-7.  Mode  Shapes  Panel  #  1  using  10  x  IS  Element  Mesh 
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It  was  nodoed  also,  from  the  convef^ence  analysis,  that  it  is  necessary  to  have  about  four 
ek^-nents  or  divisions  per  each  half-sine  wave  of  modal  di^lacement  in  order  to  d)tain  the 
precise  mode  shapes,  this  is  also  convenirat  for  the  eigenvalue  solution. 

It  must  be  pointed  out  that  there  is  an  effect  of  the  eigenvectw  on  the  natural 

frequencies  and  mode  shapes.  In  oder  to  solve  the  eigosvalue-eigaivector  proUon, 
DSHELL  uses,  as  it  was  e]q)laitwd  befrxe,  the  subq)ace  iteration  method.  This  algmithm, 
according  to  the  number  of  eigenvalues  requiered  p,  generates  the  number  of  starting 
iteraticm  vectors  q,  whoe  q>p  (actually  q  =  2p).  DSHELL  had  been  working  widi  10 
starting  iteration  vectors,  and  it  seems  that  they  were  not^  enough  for  solving  the  present 
problem,  therefore  and  adjustment  was  made  in  mder  to  obtain  good  results  and  actually 
the  number  of  starting  iteration  vecmrs  has  been  increased  to  fourteen  and  with  that 
number,  the  sdudons  obtained  were  satisfactory.  In  order  to  show  the  change  accuracy, 
the  following  eigenvectors  using  q=l0  and  14  are  conqrared  for  a  particular  shell,  as 
shown,  considering  the  higher  modes. 

One  will  notice  the  inqirovemait  for  q=  14,  furthermore,  a  table  conqrarison  of  the 
frequencies  is  also  shown.  It  becomes  apparent,  that  for  certain  cases  a  q=  14  vectors 
was  satisfactory  and  thus  this  was  used  throughout 
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( a )  Using  10  vectors 


( b )  Using  14  vectors 


Fig  3-8.  Third  mode  shiqw,  (0O/90<’]s  ply  orientation.  Panel  #  2, 
All  Edges  Sinqjly  Supported  Boundary  Condition, 
b/h  =  20  and  h/R  ==  1/20  using  «7  =  10  and  14. 


Also  die  following  taUe  shows  die  effects  of  increasing  die  number  of  vector  itnation 
fiom  <7=  10  to  9=  14  on  the  third  and  fourth  mode  for  the  especified  case. 


10  VECTORS 

14  VECTORS 

b/h  r:  10 

Mode#3 

Mode#4 

Mode#3 

Mode#4 

h/R  =  1/20 

20398.18 

20399.97 

20461.18 

20498.97 

h/R  =  1/50 

20397.15 

20398.14 

20375.15 

20419.14 

h/R  =  1/200 

20395.14 

21534.41 

20309.14 

20362.41 

1  1 

b/h  =  15 

Mode#3 

Mode  #4 

Mode#3 

Mode#4 

h/R  =  1/20 

13598.94 

13607.48 

13598.58 

13617.48 

h/R  =  1/50 

13596.58 

13598.38 

13568.24 

13586.38 

h/R  =  1/200 

13500.43 

13578.20 

13500.43 

13579.20 

1  1 

Mode#3 

Mode#4 

Mode  #3 

Mode#4 

h/R  =  l/20 

10205.80 

10297.55 

10208.94 

10397.55 

h/R  =  1/50 

9314.28 

10269.94 

9315.28 

10269.99 

h/R=  1/200 

9088.64 

10193.29 

9089.63 

10195.19 

Table  3-3.  Cmnparison  between  the  Tbird  and  Fourth  Natural  Frequencies, 
b/h  -  10,  IS  and  20,  using  ^  =  10  and  q  =  vectors. 
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IV.  Results  and  Discussions 


Hie  primaiy  olgective  of  this  thesis  was  to  exanune  the  thickness  and  curvature  effects  on 
die  (fynamic  respemse  ( natural  frequencies  and  mode  diapes  )  of  the  dicular  cylindrical, 
conqiosite  shells,  using  die  DSHELL  finite  element  program.  Also  dre  dynamic 
ciqiabiiities  of  this  code  will  be  evaluated  by  conqiaiing  to  a  reference  using  the  Galeridn 
technique. 

Previous  analyses  at  AFTT  conconing  the  dynamic  teqxmse  of  solid  cylindrical  conqxisite 
panels  have  been  developed  using  the  Galeridn  technique  and  also  including  a  higher  order 
shear  deformation  theory  ( develr^ped  Linneman  &  Pidazotto  in  1988 ). 

A  linear,  fiee  vibration  analysis  was  performed  on  the  panels  described  in  Chapter  in  in 
wdiich  danqiing  effects  were  assumed  negligiUe  and  compared  to  experiments. 

The  mode  shiqies  were  created  by  plotting  a  surface-contour  plot  of  the  eigenvector 
output  by  DSHELL.  To  help  visualize  the  mode  shapes,  orthographic  views  of  the  mode 
duqies  are  included. 

The  resuh^ng  analytical  natural  frequencies  for  this  study  and  the  ones  obtained  from  the 
Galnrkin  technique  are  compared  using  a  percent  difference  {q)proach.  In  all  cases,  the 
percoit  difference  for  a  given  mode  was  calculated  using  the  following  relation 

Pfercent  Difference  =  Wsmin  -  ^ 

foALSnON 
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where  fosiBu,  is  the  finiie  element  natund  frequency  uang  DSHELL  and  foALmn*  is  the 
natural  frequency  using  the  Galeikin  technique.  Both  frequoicies  are  given  in  Hertz  ( Hz). 

Pantl  #  1: 

This  panel  was  selected  as  a  baseline  against  which  the  model  ( using  DSHELL )  was 
cooqMired.  Since  analytical  and  experimental  results  were  available,  the  cmnpatisons  made 
he4}ed  us  to  obtained  some  important  coiM;lusions  for  further  analyses. 

The  first  consideratimi  in  analyzing  tlw  results  was  to  lodk  at  the  convergence  of  the 
eigenvalues  and  the  respective  eigenvectors.  Different  mesh  configuratirnis  were  studied: 

8  X  20, 1 1  X  12, 10  X  15,  but  only  one  ccmverged  to  the  accepted  result,  the  10  x  IS 
element  mesh  configuration  was  chosen  for  modelling  this  panel. 

Table  4-1  shows  the  predictions  provided  by  DSHELL  for  the  first  four  natural  frequencies 

using  the  36  DOF  element  and  the  10  x  IS  mesh  configuration  and  the  results  from  the 

reference  study  [7].  The  effects  of  the  degrees  of  fredom,  especially  at  the  ccxmer  nodes, 
on 

the  mode  shapes  was  also  considoed.  Hgure  4.1  shows  a  typical  comer  node  element 
used  in  this  analysis  where  the  boundary  ccmditions  are  specified. 
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Table  4-1.  Natural  Frequencies  ( Hz )  fen*  Panel  #  1 
using  the  10  x  IS  mesh  configuration 


FREE  EDGE 

32  53 

21 
20 

SIMPLE  SUPPORTED  EDGE 

Fig.  4.1  Typical  comer  tt>p  element  and  respective  node  numboing 
for  the  10  x  IS  element  mesh 


Two  sets  of  degrees  of  freedom  were  analyzed  at  each  comer  node  to  see  the  effects  of 
physically  varying  the  restraint  to  duplicate  the  actual  test  conditions.  Table  4-2  shows 
the  DOF  used  in  the  comer  element  specified  above  with  both  sets  of  DOF.  The 
behavior  of  the  oKxle  sh{q)es  were  not  affected  using  these  two  approximations. 
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B.C. 

Option 

(a) 

NODE 

NUMBER 

u 

B 

w 

HI 

HV 

m 

m 

19 

1 

1 

1 

1 

0 

1 

0 

20 

1 

1 

0 

0 

0 

0 

0 

21 

1 

1 

1 

1 

1 

1 

1 

53 

0 

0 

H 

0 

H 

0 

0 

1  n 

B.C. 

Option 

(b) 

19 

1 

1 

1 

1 

0 

1 

0 

20 

1 

1 

0 

0 

H 

H 

0 

21 

1 

1 

1 

‘  1 

0 

1 

H 

53 

0 

0 

0 

0 

0 

H 

0 

Where  the  number  1  stands  for  a  restraint  DOF 
and  0  stands  for  a  free  DOF. 


Table  4-2.  OOF  used  at  the  comer  nodes  for  Panel  #  1 


Frequency 

Number 

Experinrental 

Approach 

Mesh 

10xl5(a) 

Mesh 

10xl5(b) 

1 

102,83 

■9 

2 

mm 

199,92 

mm 

3 

209,80 

4 

IHRH 

338,05 

Table  4-3.  Effects  of  the  DOF  on  the  Natural  Frequencies  ( Hz ) 
for  options  (a)  and  (b)  shown  in  Table  4-2,  Panel  #  1 


Mod*  Nunter 

Fig.  4.2  Magnitude  of  the  Natural  Frequencies  for  the  first  four  modes, 
using  Mesh  element  10  x  IS  ( a ),  ( b )  and  reference. 


It  was  found  that  using  the  adequated  degrees  oi  freedom,  opticMi  (a)  specifically,  the 
panel  converged  in  natural  fiequencies  and  mode  shapes  to  the  result  obtained  by  using 
STAGSC-1  code  and  the  holographic  experimentation  [7].  Of  course,  absolute 
convergence  could  not  be  proven.  The  largest  percent  difference  was  7.40%  (above) 


fw  the  case  of  the  third  frequency,  and  the  smallest  was  1.68  %  ( below )  for  the  second 
natural  frequency.  The  flrst  nwde  is  followed  by  a  substantial  increase  in  frequency  to  a 
relatively  close  second  and  third  modes,  followed  by  another  substadal  increase  in 
frequency  fmr  the  fourth  mode. 

The  mode  sh{q)es  displayed  the  same  behavior  as  the  ones  shown  in  die  reference  study. 
A  graphical  representation  of  the  mode  shapes  for  this  panel  is  given  in  Figure  4.3.  By 
referring  to  this  figure,  it  can  be  seen  that  the  first  mode  was  antisymmetric,  the  second 
and  third  mode  were  symmetric,  and  the  fourth  mode  again  was  antisymmetric.  Therefore, 
trasition  from  antisymmetric  to  symmetric  modes  and  visa  versa  is  accompanied  by  a  large 
increase  in  frequency,  whereas  adjacent  modes  which  are  closely  spaced  in  frequency  were 
either  both  symmetric  or  both  antisymmetric. 
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Third  ModeShi^ 


Fourth  Mode  Shape 


Fig.  4.3  Mode  Shapes  for  Panel  #  1,  Simply  Supported-Free 
Boundary  Condition. 

It  may  also  be  observed  that  this  panel  is  deep  with  =  and  is  thin  since  S/h  = 
Therefore,  it  is  expected  that  in-plane  displacements  enter  the  eigenvector  as  important 
parameters  while  the  through  the  thickness  shear  is  unimportanL 


PanetiH: 


Using  the  DSHELL  code,  this  study  investigated  the  thickness  effects  on  the  natural 
frequency  and  mode  shapes  for  a  square  panel  varying  the  span  to  thickness  ratio  and  the 
curvature  in  order  to  measure  effects  on  two  ply  orientations,  [0®/90®]s  and  [-45®/445®]s 
in  which  each  ply  was  0.25  inches  thick,  giving  a  total  panel  thickness  of  one  inch. 

Two  boundaries  conditions  were  applied:  Simply  Supported  along  the  vertical  edges  and 
free  along  the  horizontal  edges  ( SFSF  )  and  Simply  Supported  all  around  ( SS  ). 

Different  mesh  configurations  were  also  analyze  but  the  10  x  10  element  mesh  proved  to 
be  the  most  adequate  mesh  for  the  present  study. 


NODE 

NUMBER 

u 

B 

w 

B 

m 

m 

n 

19 

1 

B 

1 

1 

0 

1 

0 

20 

1 

0 

0 

0 

0 

0 

0 

21 

1 

1 

1 

1 

1 

1 

1 

0 

1 

0 

0 

0 

0 

0 

'HHHi 

0 

1 

1 

0 

1 

0 

1 

Table  4-4.  Typical  comer  top  element  and  the  characteristic  node  numbering 

for  the  10  x  10  mesh  configuration 
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Table  4-4  shows  the  DOF  used  in  a  typical  comer  element  fw  this  panel,  and  Figure  4.4 
shows  the  specified  element  (comer  top)  above  and  the  respective  nodes  numbering. 
The  nodes  with  astmks  in  that  figure  are  only  considered  in  the  SS  boundary  condition. 

FREE  EDGE 

32  53 


SIMPLE  SUPPORTED  EDGE 

Fig.  4.4  Typical  top  comer  element  and  the  respective  node  numbering 

for  the  10  x  10  mesh  ccmfiguration 

The  panels  studied  were  square,  a  =  b  =  1,  where  a  and  b  are  the  edges  along  the  x  and 
s  axis  respectively;  the  thickness  was  assumed  ccmstant  and  for  the  boundary  conditions 
specified,  different  ratios  of  thickness  to  radius  of  curvature,  h/R,  were  assumed  ( 1/20, 
l/SO  and  1/200 )  and  specific  values  of  circumferential  length  to  thickness  ratios,  b/h  ( 10, 
20, 30, 40  and  50 )  were  considered. 
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AU  E^es  SimfAy  Supported  Boundary  Condition, 

Figure  4.5  shows  the  effect  of  curvature  h/R  on  the  natural  frequencies  for  panel  #  2, 
[(FV90O]s  ply  (Hientation,  with  simply  supported  boundary  condition. 


(b/hx10) 

Fig  4.5  Curvature  Effects  on  the  First  Natural  Frequency,  Panel  #  2, 

Simply  Supported  Boundary  Condition,  [0'^/90^]s  laminate 

Some  important  trends  are  identified;  the  natural  frequencies  are  seen  to  increase  as  h/R 
is  increased  (smaller  radius  of  curvature).  Also  the  effect  of  increasing  the  span  to 
thickness  ratio,  b/h,  (making  the  panel  deeper)  is  seen  to  lower  the  frequencies.  In 
addition,  for  relatively  "close"  Boundary  Conditions  (shallow  panels),  b/h  =  10,  a  small 
difference  in  the  magnitude  of  the  first  frequency  is  shown  for  the  three  span  to  thickness 


ratios  studied;  this  behavior  changes  for  relatively  "further  apart"  boundaries,  b/h  >  40, 
where  the  large  differences  in  frequency  appeared.  This  implies  that  shallow  panel 
response  is  insensitive  to  radius  of  curvature  but  as  the  panel  becomes  deeper  the 
frequency  is  affected  the  curvature. 

Hgure  4.6  shows  the  effects  of  the  radius  of  curvature  on  the  natural  frequency  for  the 
[-45^/-f450]s  laminate.  The  fundamental  frequencies  are  seen  to  increase  as  h/R  is 
increased  but  die  frequencies  decrease  by  increasing  the  span  to  thickness  ratio. 


Fig  4.6  Curvature  Effects  on  the  First  Natural  Frequency,  Panel  #2, 
Simply  Supported  Boundary  Condition,  [-45^/+  45^]$  laminate 
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Hie  small  difference  in  frequencies  for  "close”  boundaries  as  seen  in  the  last  case,  Hg  4.5, 
is  present  fw:  the  [-45®/+4501s  ply  (mentation  for  h/R  >  1/50  ( Fig  4.6 ).  The  larger  h/R 
responded  with  the  largest  frequency  no  matter  what  shallowness  parameter  tVh  was 
considered  From  the  figures  above  can  be  seen  also  that  for  shallow  panels  (6/c  very 
small),  the  S/h  implies  thick  panel  and  the  effects  of  throu^  the  thickness  shear  is  more 
noticed. 

Cooqiarisons  can  also  be  made  between  the  natural  fiequencies  obtained  for  the  two 
laminates.  Hgure  4.7  shows  the  tendency  of  the  curves  for  the  fundamental  frequency 
fiv  die  laminate  orientations  described  the  value  h/R  =  20  wa^  chosen  for  this  analysis. 


Fig.  4.7  CcMnparison  for  the  Fundamental  frequency  between  the  [0^/90^]$  and 

[-45®/+45®]s  laminates.  Simply  Supported  Boundary  Condidonh/R  =  1/20. 
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Both  curves  show  the  same  behavior,  the  magnitude  the  first  fitequracy  is  higher  for  the 
[-4S^/-f450]s  than  the  [0°/90‘’]s  laminate  since  that  ply  (xientation  makes  the  panel  stiffer. 
This  is  because  of  the  frequency  equation  depends  on  the  extensional  ( )  and  holding 
coefficients  (  Dy  ),  [18:344],  and  they  axe  higher  for  the  [-45®/+45®]s  than  the  10®/90®]s 
laminate.  Also,  by  increasing  the  span  to  thickness  ratio,  the  first  frequency  decreases  in 
bodi  cases.  In  addition,  for  (b/h)  <  30,  the  variation  between  the  two  frequencies  is 
significant,  but  as  b/h  increases  the  conqiarison  frnr  the  two  lay-iqis  convoge. 

Up  to  this  point  these  results  are  not  new.  Unneman  [8]  examined  symmetrical  laminates 
using  two  theories:  Galeridn  Technique  aikl  a  Higher  Order  Shear  Theory;  the  validation 
of  the  results  in  this  analysis  was  obtained  frtxn  comparistm  to  his  results. 

The  values  fOT  the  first  natural  frequency  uring  DSHELL  and  the  Galeridn  technique  for 
the  ratio  h/R  =  20  are  shown  below.  The  difference  percent  is  also  introduced.  Table  4-5 
shows  the  first  natural  frequency  far  this  panel,  [0^/90^]s  laminate,  with  simply  supported 
boundary  condition.  The  highest  difference  obtained  between  the  Galeridn  technique  and 
the  DSHELL  approximations  was  1 1.34  %  fra*  the  ratios  b/h  =  30  and  h/R  =  1/50. 

TaUe  4-6  also  shows  the  first  natural  frequency  fOT  the  [-45^/445^]^  laminate,  same 
boundary  condition;  the  highest  difference  obtained  between  the  Galeridn  technique  and 
the  DSHELL  approximations  was  9.40  %  for  the  ratios  b/h  =  50  and  h/R  =  1/20,  while 
the  lowest  percent  in  difference  was  1 .08  %.  The  percent  differences  are  relatively  small, 
thoefOTe  DSHELL  with  the  associated  model  was  used  for  subsequent  conditions. 
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ll/hr:  10 

HSSSfl 

DSHELL 

FEMcode 

Iterance 

Percent 

h/R»l/20 

10800 

10702.7 

+  0,90% 

h/R=l/S0 

IQSOO 

10482,2 

+  0.17  % 

h/R»  1/200 

NA 

10441,4 

- 

II 

Gakridn 

TechniQue 

DSHELL 

FEMcode 

Difiference 

Percent 

h/R=l/20 

4200 

4016,3 

+  4,36% 

h/R=l/50 

3200 

3124.4 

+  2,36% 

h/R=  1/200 

N.A 

2919,2 

- 

li/h  =  30 

DSHELL 

FEMcode 

Difference 

Percent 

h/R=l/20 

2950 

+  6,43% 

h/R=l/50 

hei^shh 

1773,2 

+  11,34% 

h/R«  1/200 

NA 

1355,5 

• 

h/hs40 

DSHELL 

FEMcode 

Difference 

Percent 

h/R=l/20 

1900,3 

+  4,98% 

h/R=l/50 

1450 

1400,2 

+  3.43  % 

h/R=  1/200 

• 

809,3 

- 

b/h  =  50 

Galeridn 

Technique 

DSHELL 

FEMcode 

Difference 

Percent 

h/R=l/20 

1700 

1575,8 

+  7,30% 

h/R=l/50 

1100 

1071,2 

+  2,62  % 

h/R=  1/200 

N.A 

570,7 

- 

(N.  A  =  Not  AvalaUe) 

Table  4-5  Comparison  betwera  l^HELL  and  the  Galeridn  Technique  for 
the  Fundamental  Frequency,  Panel  #  2,  [0°/90^]s  ply  orientation 
Simply  Supported  Boundary  Ccmdition. 
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DSHELL 

FEMcode 

Diffeience 

Percent 

h/R  =  l/20 

14300 

13672,2 

+  4,39% 

h/R=l/50 

12900 

12349,6 

+  4,26% 

h/R^^  1/200 

m 

12092.0 

- 

h/hslO 

Qakridn 

Technique 

SHELL 

FEMcode 

Differmce 

Percent 

h/R  =  l/20 

7700 

7285,9 

+  5,38  % 

h/R=l/50 

4500 

4389,5 

+  2,46% 

h/R=  1/200 

N.A 

3504,1 

- 

b/hs30 

DSHELL 

FEMcode 

Diffeience 

Percent 

h/R=l/20 

3668,9 

+  8.28  % 

h/R=l/S0 

3200 

3165,3 

+  1,08% 

h/R=  1/200 

N.A 

1720,8 

• 

DSHELL 

FEMcode 

Differmce 

Pucent 

h/R  =  l/20 

2550 

2414,9 

+  5,30% 

h/R=l/50 

2150 

2107,7 

+  1.97% 

h/R  =  1/200 

N.A 

1135,6 

- 

b/h  =  50 

Galerkin 

Technique 

DSHELL 

FEMcode 

Diffnence 

Percent 

h/R=l/20 

2100 

1902,7 

+  9,40% 

h/R  =  l/50 

1500 

1463,6 

+  2,43% 

h/R=  1/200 

N.A 

906,1 

- 

(NA  =NotAvailibk) 

Table  4-6  Comparison  between  DSHELL  and  the  Galeridn  Technique 

for  the  First  Frequency,  Panel  #  2,  [-45®/445®]5  ply  orientation 
Simply  Supported  Boundary  Condition. 


4-15 


Although  only  die  fundamental  fitequency  has  been  included  in  the  previous  analysis,  the 
eigenvectors  cone^cniding  to  die  first  four  natural  fiequencies  for  the  [0^/90^]$  and 
[-4S^/+4SO]s  laminate,  for  b/h  =  20,  with  Sinqily  Supported  all  edges  Boundary  Gxidition 
are  included  in  the  next  Hgures.  As  stated  in  the  modelling  chtqiter,  a  subqiace  of  14 
vectors  was  used  throughout  with  a  cmivergence  tolerance  on  eigenvalue  of  lO'*^ ,  thus 
twelve  digit  accuracy  and  six  digit  accuracy  for  the  vector  conqionents  of  the  eigenvector 

-  A,*| 

using  a  Raylei^  quotient  relation  [21].  The  c(»veigence  equation  is  '  *  ^  tol 

where  k  is  the  given  iteration. 

Hg  4.8  shows  the  mode  shqie  for  the  [0°/90^]s  laminaiB,  under  die  sinqily  supported  all 
edgesboundaryconditi(Mi,forthespecific  values  of  h/R=:  1/20  and  b/h  =  20.  The  first 
mode  is  symmetric,  the  second,  dmd  and  fourth  mode  are  antisymmetric. 
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Third  Mode  Shape  Fourth  Mode  Shape 

Fig.  4.8  Mode  Shapes  for  Panel  #  2,  Singly  Supported  Boundary 
Condition,  [O^Wls h/R  =  1/20,  b/h  =  20 

Bg  4.9  shows  the  mode  shape  for  the  [■450/+450]s  laminate,  under  the  simply  supported 
boandaiycondition,forthespecific  values  of  b/R=  1/20  and  h/h  =  20.  The  first  mode  is 
symmetric;  the  second  and  diird  and  fourth  inode  ate  antisymmetric. 
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Rrst  Mode  Shape 


Second  Mode  Shq)e 


Fig.  4.9  Mode  Slu4>es  for  Panel  #  2,  Smq>ly  Supported  Boundary 
G>ndition  [-45®/450]s  landnate,  h/R  =  1/20,  b/h  =  20 
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A  complete  set  d  results  ai  this  analysis  incluiding  the  first  four  natural  fitequoicies,  fw 
both  ply  orientation  and  the  specified  boundary  conditions,  is  given  in  the  following 
tables. 


h/hiilO 

Mode  #1 

Mode  #2 

Mode  #3 

Mode  #4 

h/R=l/20 

10702.72 

18977.13 

20461.18 

20498.97 

h/R=l/50 

10482.24 

18649.47 

20375.15 

20419.14 

h/R=  1/200 

10441.45 

18542.59 

20309.14 

20362.41 

1  1 

h/hsl5 

Mode  #1 

Mode  #2 

Mode  #3 

Mode#4 

h/R=l/20 

5685.29 

9374.52 

13598.58 

13617.48 

h/R=l/50 

5127.56 

9360.28 

13568.24 

13586.38 

h/R=  1/200 

5016.13 

9266.70 

13500.43 

13549.20 

1  1 

1  b/hs20 

Mode  #1 

Mode  #2 

Mode  #3 

Mode  #4 

s 

II 

4016.82 

5524.13 

10205.80 

10397.55 

h/R=:l/50 

3124.43 

5497.42 

9314.28 

10269.94 

h/R»  1/200 

2919.27 

5429.79 

9088.64 

10193.29 

1  1 

b/hr:30 

Mode  #1 

Mode  #2 

Mode  #3 

Mode  #4 

h/R=l/20 

2760.3 

3036.1 

4904.9 

5698.4 

h/R=l/50 

1773.2 

2563.2 

4754.5 

4988.8 

h/R=  1/200 

1355.5 

2524.7 

4417.8 

4708.8 

1 

b/h  =  40 

Mode  #1 

Mode  #2 

Mode  #3 

Mode  #4 

h/R=l/20 

1900.36 

2702.33 

2842.59 

4030.63 

h/R=l/50 

1400.27 

1524.24 

2737.03 

3137.14 

h/R=  1/200 

809.38 

1442.29 

2613.72 

2836.48 

r 

II 

1 

Mode  #1 

Mode  #2 

Mode  #3 

Mode  #4 

h/R=l/20 

1575.89 

1931.45 

2477.51 

3127.58 

h/R=l/50 

1071.26 

1266.30 

1853.52 

2233.47 

h/R=  1/200 

570.72 

935.18 

1757.98 

1838.85 

Table  4-7  First  four  Natural  Frequencies,  Panel  #2,  Simply 

Supported  Boundary  Condition,  [00/90®]  g  ply  orientation. 
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h/h-lO 

Mode  #1 

Mode  #2 

Mode  #3 

Mode  #4 

h/R=l/20 

13672.26 

22404.69 

26239.22 

29104.94 

h/R=l/50 

12349.65 

22138.28 

25899.37 

28886.96 

h/R=  1/200 

12092.08 

22076.65 

25819.76 

28839.93 

b/h=lS 

Mode  #1 

Mode  #2 

Mode  #3 

Mode  #4 

h/R=l/20 

8840.34 

11772.55 

14737.46 

18408.56 

h/R  =  1/50 

6468.87 

11367.77 

13818.49 

17716.12 

h/R=  1/200 

5919.98 

11232.03 

13698.42 

17578.29 

1 

Mode  #1 

Mode  #2 

Mode  #3 

Mode  #4 

h/R  =1/20 

8285.94 

7427.99 

10384.23 

12037.38 

h/R  =1/50 

4389.53 

6935.45 

8615.55 

11019.53 

h/R  =  1/200 

3504.12 

6717.43 

8379.37 

10777.59 

Mode  #1 

Mode  #2 

Mode  #3 

Mode  #4 

h/R  =  1/20 

3668.95 

6254.74 

6431.54 

8701.83 

h/R  =1/50 

3165.36 

3459.07 

4667.79 

5666.27 

h/R  =1/200 

1720.89 

3167.73 

4026.06 

5187.10 

Mode  #1 

Mode  #2 

Mode  #3 

Mode  #4 

h/R  =1/20 

2414.96 

3685.05 

5469.36 

5625.56 

h/R  =1/50 

2107.75 

2872.94 

3512.95 

3618.75 

h/R  =1/200 

1135.68 

1859.18 

2373.24 

3051.14 

b/h  =  SO  Mode  #1  Mode  #2  Mode  #3  Mode  #4 
h/R=l/20  1902.79  2678.49  3873.12  4603.97 

h/R  =  l/50  1463.64  2356.94  2756.15  3119.88 

h/R=  1/200  906.12  1245.63  1605.09  2039.39 


Table  4-8  First  four  Natural  Frequencies,  Panel  #  2,  Simply  Supported 
Boundary  Condition,  [-45®/+  45®]s  ply  orientation. 


The  next  figures  show  the  behavior  of  the  second,  third  and  fourth  mode.  They 
constructed  using  the  data  above. 
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Hgs  4.10  and  4. 1 1  show  the  behavior  of  the  frequencies  for  the  second  mode  under  the 
boundary  condition  considered.  According  to  these  figures,  the  frequencies  for  the 
[-45®/445®]s  laminate  are  higher  than  the  ones  for  the  [0®/90®]s  ply  orientation.  The 
panel  stiffness  is  higher  in  the  [-4SO/+45^]s  laminate  as  indicated  previously.  Also  the 
effect  of  the  through  the  thickness  shear  becomes  evidently  higher  for  the  [-45®/4450]5 
laminate. 

Also  from  fig  4.1 1,  for  values  b/h  >  20  (deeper  and  thinner  shells)  through  the  thickness 
shear  effects  decreased  but  the  curvature  effects  becomes  significant.  Also  for  deeper 
shells,  b/h  >  50,  on  the  same  value  of  curvature,  the  behavior  of  the  curve  seems  to  reach  a 
specific  value,  for  both  orientations,  this  value  should  correspond  to  the  frequency  of  a 
shell  with  no  shear. 


Fig  4. 10  Curvature  Effects  for  the  Second  Natural  Frequency,  Panel  #  2, 
Simply  Supported  Boundary  Condition,  [0®/90O]s  laminate. 
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Fig  4. 11  Qirvature  Effects  for  the  Second  Natm^  frequency.  Panel  #  2 
Simply  Suppmtol  Boundary  Condition,  [-45®/+4501s  laminate 


Figs  4.12  and  4.13  show  the  behavior  of  the  frequencies  for  the  third  mode.  According  to 
these  figures  the  frequencies  for  the  [-4SO/44S^]5  laminate  are  higher  than  the  ones  for 
the  [0®/90®]s  ply  orientation.  The  effect  of  the  through  the  thickness  shear  is  also 
evidently  higher  in  the  [-45®/+45®js  ®  10®/90®]s  laminate.  The  effect  of 

curvature  is  practically  negligible  for  the  [0®/90®]s  laminate  ( since  the  curves  are  very 
close,  but  in  the  [-45®/445®]s  laminate,  it  becomes  important 

The  same  trends  are  also  seen  in  figs.  4.14  and  4.15,  for  the  fourth  natural  frequency. 
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Hg  4. 12  Curvature  Effects  for  the  Third  Natural  Rnequency,  Panel  #  2, 
Simply  Supported  Boundary  Condition,  [(P/9(P]s  laminate 


Fig  4.13  Curvature  Effects  for  the  Third  Natural  frequency,  Panel  #2 
Simply  Supported  Boundary  Condition,  [-45®/445®]s  laminate 
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0  1  2  3  4  5  6 

b/h 


Fig  4.14  Curvature  Effects  f(»^  the  Fourth  Natural  Frequency,  Panel  #2 
Simply  Supported  Boundary  Condition,  [00/90^]$  laminate 
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Hg4.15  Curvature  Effects  for  the  Fourth  Natural  Frequency,  Panel  #2 
Simply  Supported  Boundary  Condition,  [-45®/+45®]s  laminate 
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Also  comparisons  can  be  made  between  the  ply  orientatimis  fm-  the  specific  value  of 
h/R  =  1/20,  Figs  4.16, 4.17  and  4.18  show  the  behavimr  of  the  second,  third  and  fourth 
mode  respectively. 

It  is  noticed  that  the  frequency  magnitude  for  the  [-4S^/44SO]s  ply  (mentation  is  higher 
than  the  corresponding  values  for  the  [(P/90O]s  laminate  for  the  second,  third  and  fourth 
noodes.  In  addition,  for  b/h  >  30,  the  difference  in  frequency  between  both  orientations 
is  veiy  small. 


Fig  4.16  Qmiparisons  between  [0O/90®]s  and  [-45®/+45®]s  laminate.  Panel  #  2 
Second  Natural  Frequency,  Simply  Supported  Boundary  Condition 
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Fig  4.17  Ounparisons  between  [(P/90°]s  and  [-45®/+45®]ff  laminate.  Panel  #  2 
Third  Natural  Frequency,  Simply  Supported  Boundary  Condition 


Fig  4. 18  Comparisons  between  [0O/90®]s  and  [-45®/+45®]s  laminate.  Panel  #  2 
Fourth  Natural  Frequency,  Simply  Supported  Boun-^  Condition 
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Simpfy  Supported-Free  Boundary  Condition. 


The  next  case  investigated  was  the  Singly  Suppcmed-Free-Simply  Supported-Free 
Boundaiy  Condition.  The  DOF  applied  in  this  boundary  condition  were  q)ecified  in 
Table  4-4.  DSHEUL  again  showed  flexibility  in  getting  the  eigenvalues  and  the 
cmresponding  noode  shapes  for  this  boundary  oHidition,  so  the  calculations  were  straight 
forward. 

Following  the  pattern  used  with  the  previous  boundary  condition,  the  same  thickness  to 
radius  of  curvature  ratios,  h/R,  were  used  (1/20, 1/50, 1/200)  for  specific  values  of 
circumferential  length  to  thickness  ratios  b/h  (10, 20, 30, 40  and  50 ). 

The  next  taUes,  Table  4-9  and  Table  4-10  respectively,  show  the  magnitude  of  the  first 
four  natural  frequencies  using  DSHELL  finite  element  program  for  both  ply  orientation 
t0®/90Ols,and  [-450/4450]s. 
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UiBliMUl 

Mode  #1 

Mode  #2 

Mode  #3 

Mode  #4 

s 

II 

9043.4 

15759.8 

16582.6 

h/R=l/50 

5240.8 

7259.0 

15240.6 

16386.6 

h/R=  1/200 

5063.2 

6875.7 

15226.7 

16334.4 

1  1 

Mode  #1 

Mode  #2 

Mode#3 

Mode  #4 

h/R=l/20 

3961.8 

6688.1 

8052.4 

9365.8 

h/R=l/50 

2764.7 

4129.6 

7896.7 

9256.8 

h/R=  1/200 

2408.6 

3393.3 

7605.5 

9152.6 

1  1 

Mode  #1 

Mode  #3 

Mode  #4 

h/R=l/20 

3181.3 

4469.7 

5575.2 

57%.2 

h/R=l/50 

1941.3 

3125.6 

4862.5 

5504.1 

h/R=  1/200 

1431.2 

2075.9 

4745.3 

5490.3 

1  1 

ISiSlilUl 

Mode  «1 

hfode  #2 

Mode#3 

Mode  #4 

h/R=l/20 

1965.3 

2497.1 

2964.8 

3910.8 

h/R=l/50 

1375.8 

2162.0 

2494.3 

2676.8 

h/R=  1/200 

731.4 

1122.4 

2197.5 

2612.5 

1  1 

Mode  #1 

Mode  #2 

Mode  #3 

Mode  #4 

h/R  =  l/20 

1162.4 

2033.7 

2104.3 

2585.0 

h/R=l/50 

1054.2 

1228.2 

1654.0 

2162.4 

h/R=  1/200 

502.5 

819.1 

1259.1 

1530.1 

1  1 

lESSil^ 

Mode  #1 

Mode  #2 

Mode  #3 

Mode  #4 

h/R=l/20 

787.8 

1503.3 

1688.0 

1720.2 

h/R=l/50 

636.4 

1047.4 

1327.9 

1527.7 

h/R=  1/200 

404.2 

699.0 

814.7 

1011.8 

Table  4-9  First  Four  Natural  Frequencies  using  DSHELL,  Panel  #  2, 

[0^/90^]s  laminate.  Singly  Supported-Free  Boundary  Condition. 


As  expected,  by  using  this  boundary  condition,  the  frequencies  obtained  ate  lower  than  the 
values  for  the  previous  boundary.  For  bodi  laminates ,  [0^/90*’]$  and  [-45®/+45®]s, 
die  values  of  fiequencies  obtained  are  about  50%  the  corre^nding  value  of  frequency 
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calculated  for  the  Simply  SuiqxMted  all  edges  boundary  cmdition. 


lEmnoi 

Mode  #1 

Mode  #2 

i7>7y?yf'i 

Mode  #4 

h/R  =>1/20 

6851.8 

10775.6 

18488.0 

18814.4 

h/R  =  1/50 

5865.1 

9824.4 

18172.9 

18412.6 

h/R  =1/200 

5665.9 

9644.1 

18127.6 

18287.9 

1  1 

Mode  #1 

Mode  #2 

Mode  #3 

Mode  #4 

h/R  =1/20 

4786.3 

6671.5 

9457.4 

10410.5 

h/R  =1/50 

3205.8 

5110.7 

9250.4 

9352.1 

h/R  =1/200 

2805.9 

4741.0 

9111.4 

9304.2 

1 

13iSl>X4)l 

Mode  #1 

Mode  #2 

Mode  #3 

Mode  #4 

h/R  =  1/20 

4147.1 

4836.3 

6145.7 

7085.6 

h/R  =1/50 

2306.3 

3415.1 

5671.7 

5852.6 

h/R=  1/200 

1702.4 

2835.8 

5455.7 

5690.7 

1  1 

ESTBIXGl 

Mode  #1 

Mode  #3 

Mode  #4 

h/R=l/20 

2310.2 

3411.5 

4290.7 

4353.3 

h/R  =1/50 

1767.3 

2224.6 

2865.7 

3263.9 

h/R=  1/200 

883.8 

1402.5 

2616.1 

2734.4 

1  1 

isisiiXyi 

Mode  #1 

Mode  #2 

Mode  #3 

Mode  #4 

h/R  =1/20 

1448.2 

2526.4 

2916.6 

3175.3 

h/R  =  1/50 

1262.2 

1586.9 

2194.4 

2232.0 

h/R  =  1/200 

612.5 

906.2 

1582.9 

1608.2 

1  1 

b/h  =  SO  (6) 

Mode  #1 

Mode  #2 

Mode  #3 

Mode  #4 

h/R  =1/20 

951.5 

1757.5 

1990.1 

2413.3 

h/R  =  1/50 

760.8 

1427.0 

1713.0 

1828.2 

h/R=  1/200 

503.9 

692.8 

1067.5 

1112.9 

Table  4-10  First  Four  Natural  Frequencies  using  DSHELL,  Panel  #2, 

[-45®/445®]s  laminate.  Simply  SuppOTted-Free  Boundary  Condition 


The  fundamental  mode  will  be  discussed  next  since  it  is  most  inqxntant  in  the  analysis. 
Hguie  4.19  shows  the  effect  of  curvature  h/R  on  the  fundamental  frequency  for  panel  #  2, 
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[0^/90O]g  ply  orientati<Mi.  As  it  is  seen,  as  the  curvature  increases,  the  panel  bec<Hnes 
deeper  and  sdffer  and  the  natural  hequencies  increase. 


Hg  4. 19  Curvature  Effects  on  the  Fundamental  Frequency,  [0^/90°]$ 
ply  (mentation.  Simply  Supported-Free  Boundary  Condition 


Some  tendencies  are  also  recognized.  The  natural  fiequencies  are  seen  to  increase  as  h/R 
is  increased  due  to  the  mranbrane  and  bending  coupling.  Also,  the  effect  of  increasing 
the  span  to  thickness  ratio,  bi/h,  is  seen  to  lower  the  fiequencies.  The  same  behavior  was 
noticed  ftff  the  [-45®/+45®]s  laminate,  Bguie  4-20.  Both  figures  present  basically  the 
same  tendencies  except  that,  the  curves  for  the  [-45^/445^]^  laminate  are  higher  in 
magnitude  than  the  corresponding  curve  for  foe  [00/90^]$  laminate  as  indicated 
previously. 
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Rg  4.20  Curvature  Effects  on  the  Natural  Frequencies,  Panel  #  2,  [-45®/445®]5 
ply  orientaticMi,  Simply  Suj^rted-Free  Boundary  Condition 


The  next  figure  shows  the  difference  in  natural  frequencies  for  both  laminates.  Again,  the 
natural  frequencies  for  the  laminate  are  higher  than  the  value  for  the 

[0P/90^]s  orientation  due  to  the  higher  stiffness  and  thickness  effects  in  the  first  layup. 
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Fig.  4.21  Comparison  in  the  First  Natural  Frequency  between  the 
[0®/90®]s  and  [-450/+45®]s  ply  orientations,  Sinq)ly 
Supported-Free  Boundary  Condition,  h/R  =  1/20 


Also  from  this  figure  it  can  be  seen  that  for  values  b/h  <  30,  the  variation  between  the  two 
fiequencies  is  significant  but  as  b/h  increases,  the  comparison  for  both  laminates  convo-ge. 

Aldiough  (xily  die  first  natural  finequency  has  been  included  in  the  previous  analysis  (fig 
4.19),  the  mode  shapes  corresponding  to  the  first  four  natural  frequency  for  the  [0^/90O]s 
and  [-4SO/+45^]s  ply  (mentation  with  simply  supported-free  boundary  condition  are 
included  in  the  next  figures. 
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Frcnn  figure  4.22,  it  can  be  seen  that  the  first  nxxie  sh^  for  the  [0'^/90O]s  ply  orimtation, 
is  symmetric,  the  second,  third  and  fourth  mode  are  antisymmetric. 


First  Mode  Shape  Second  Mode  Shape 


Third  Mode  Shtq)e  Fourth  Mode  Shape 


Hg  4.22  Mode  Shtqpes  fw  Panel  #  2,  [0®/90O]s  ply  orientation,  b/h  =  20, 
h/R  =  1/20,  Simply  Supported-Free  Boiuidary  Condititn. 
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From  figure  4.23,  it  can  be  seen  that  the  first  mode  shi^  for  the  [-4SO/-i45°]g  ply 
orientation,  is  symmetric,  dw  second,  diird  and  fourth  mode  are  antisymmetric,  is  the  same 
as  the  [0^/90°]$  ply  orientation. 


Third  Mode  Shape  Fourth  Mode  Shape 


Rg  4.23  Mode  Shapes  fiir  Panel  #  2,  [-450/+450]  ply  orientation,  b/h  =  20, 
h/R  =1/20,  Simply  Supported-Ftee  BoundaryCondition 


In  the  following  figures  the  curvature  effects  on  the  second,  third  and  fourth  ftequencies 
for  both  ply  orientation  are  shown.  Figures  4.24  and  4.25  present  the  curvature  effects 
on  the  second  natural  frequency  for  both  laminates  respectively.  For  both  laminates,  the 
frequency  is  seen  to  increase  as  h/R  increases,  but  the  frequency  decreases  by  increasing 
the  distance  between  the  boundaries.  Also  the  frequencies  for  the  [-45®/+45®]s  ply 
orientation  are  higher  than  the  corresponding  values  for  the  [0^/90°]  g  laminate,  for  the 
four  frequencies,  through  all  interval  b/h  analyzed. 


Fig  4.24  Curvature  Effects  on  the  Second  Natural  Frequency,  [0®/90®]s 
ply  orientation.  Simply  Supported-Free  Boundary  Condition 
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In  addition,  the  difference  in  frequency  among  the  three  ratios  of  curvature  seems  to 
remain  constant  for  each  ply  orientation  along  the  interval  b//h  study,  although  fw  closed 
boundaries,  b/h  <  20,  this  difference  is  larger  in  the  [0^/90°]5  ply  orientation.  If  these 
figures  are  compared  with  the  respective  ones  obtained  for  the  previous  boundary 
( figures  4. 10  and  4. 1 1 ),  it  is  seen  that  the  difference  in  frequency  among  the  three  ratios 
h/R  studied,  for  each  b/h  value,  is  higher  in  the  present  boundary  condition  case  for  both 
laminates  than  the  corresponding  ply  (mentation  in  the  previous  boundary. 


Fig  4.25  Curvature  Effects  on  the  Second  Natural  Frequency,  [-45®/+45®]s 
ply  orientation.  Simply  Supported-Iree  Boundary  C(Hidition 
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Hgures  4.26  ami  4.27  show  the  curvature  effect  on  the  third  natural  frequency  for  both 
laminates  respectively;  frequency  increases  as  h/R  increases  too,  but  it  decreases  by 
increasing  the  distance  between  boundaries.  Also,  the  effect  of  the  boundary  condition  is 
noticed,  by  comparing  these  figures  with  figures  4.12  and  4.13,  corresponding  to  the  same 
frequency  but  simply  supported  all  edges,  the  frequency  magnitudes  in  the  present  case  are 
lower  than  the  values  for  frequencies  obtained  for  the  previous  boundary. 


Fig  4.26  Curvature  Effects  on  the  Third  Natural  Frequency,  [0°/90®]s 
laminate.  Simply  Supported-Free  Boundary  Condition 
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Fig  4.27  Curvature  Effects  on  the  Third  Natural  frequency,  [-45®/+45®]s 
ply  orientation.  Simply  Supported-Free  Boumlary  Condition 


Hgures  4.28  and  4.29  show  the  curvature  effect  tm  the  fourth  natural  frequency  for  both 
laminates  re^)ectively;  frequency  increases  as  h/R  increases  too,  but  it  decreases  by 
increasing  the  distance  between  boundaries.  These  curves  follow  exactly  the  same  trend 
seen  in  the  third  mode,  the  effects  of  curvature  are  not  observable  at  the  shallower  end  of 
the  composite  (b/h  <  20),  specially  in  the  [0^/90O]s  ply  orientation,  but  diey  are  more 
evident  at  the  deeper  end  in  the  [-45®/445®]s  ply  (mentation. 
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Fig  4.28  Curvature  Effects  on  the  Fourth  Natural  Frequency,  [0^/90^]$ 
laminate.  Simply  Supported-Free  Boundary  Condition 


Fig  4.29  Curvature  Effects  on  the  Fourth  Natural  Frequency,  [-450/445®]  g 
ply  orientation.  Simply  Suppcvted-Free  Boundary  Condition 
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Following  the  anafysis  made  with  the  previous  boundary  condition,  comparisons  between 
both  ply  orientations  can  also  be  made  for  h/R  =  1/20. 

Hgures  4.30, 4.31  and  4.32  show  the  difference  between  the  [0®/90Ols  and  [-45®/4450]s 
ply  (mentation  fen*  the  second,  third  and  fourth  natural  frequency  for  the  simply  supported- 
fiee  boundary  condition.  It  is  seen  that  the  natural  frequencies  for  the  [-45°/-i450]s 
laminate  are  higher  than  the  corresponding  frequencies  for  the  [0^/90°]$  ply  (mentation, 
this  difference  is  about  20%,  for  the  three  cases,  and  it  remains  approximately  constant 
through  all  the  interval  b/h  analyzed.  Also,  it  is  evident  that  the  effect  (tf  the  through  the 
thickness  shear  on  the  [-450/+4S°]s  ply  orienttuion  is  greater,  spettially  at  the  shallower 
end  of  die  composite,  although  this  effect  bectxnes  negligible  as  the  deepness  is  increased 


Hg  4.30  Comparison  between  [0®/90®]s  and  [-45®/+450]s  ply  orientation.  Second 
Natural  Bequency,  Simply  Supported-Free  Boundary  Condition 
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Hg4.31  G)nq)arison  between  [00/90®]  s  and  [-45®/+45®]s  laminate.  Third 
Natural  Frequency,  Sin^)l}  •'  *pported-Free  Boundary  Condition 


Fig  4.32  Om^Muison  bciweoi  [0®/90®]5  and  [-450/+45®]s  ply  orientation. 

Fourth  Natural  Frequency,  Simply  Supported-Free  Boundary  Condition 
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V.  ConduskMU  and  Recomendatioiis 


An  analytical  study  was  carried  out  using  the  DSHELL  finite  element  program  in  the 
solution  of  dynamic  response,  natural  frequencies  and  mode  shapes,  for  curved  Gr/Ep 
con^x)site  panels  considering  through  the  thickness  shear  strain. 

Based  upon  on  the  analysis  developed,  the  fr^owing  conclusions  were  reached. 
Con^Mtisons  for  a  panel,  used  as  a  baseline  study,  with  previous  experimentation  using 
Hdognqthic  Interferometry,  and  analytical  Judies  using  STAGSC-1  finite  element  code, 
proved  to  be  excellent  The  frrst  four  natural  frequencies  and  the  respective  mode 
shapes  obtained  in  this  study  matched  completely  with  the  results  in  the  previous  analysis. 

It  was  found  that  modelling  the  conqtosite  shell  panel  was  dictated  by  improving  the 
accuracy  of  natural  frequencies  and  mode  sh{q)es.  The  result  from  this  was  that  four 
nodes  have  to  be  positioned  per  each  half  sine  wave  framed  and  the  size  of  the  subspace 
in  the  eigenvalue  solution  was  fourteen  vectors  in  arriving  at  the  four  eigenvalues 
reqi^ted. 

Good  results  were  obtained  for  a  given  panel,  oxisidering  the  two  boundary  conditions 
studied,  all  edges  simply  supported  and  simply  suppmted-free.  The  output  obtained 
fixan  DSHELL  was  close  to  the  results  obtained  fnnn  the  Galeridn  Technique.  The 
accuracy  and  flexibility  of  the  program  was  proved. 

The  curvature  of  the  cmnposite  shells  were  varied  as  was  the  ^>an  to  thickness  ratio  in 
order  to  measure  the  effects  for  each  ply  orientation,  under  the  specified  boundary 
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conditioiis.  The  remits  showed  that,  as  the  shell  becomes  deepor  the  frequency  becomes 
smaller. 

Hndings  show  that  larger  curvature  increases  the  natural  frequencies  for  both  laminates 
[0P/90O]^  and  [-4S^/44S*’]s  due  to  the  fact  diat  the  in-plane  displacements  (u  and  v 
q)ecially)  becomes  more  important  as  the  curvature  increases.  This  was  found 
considering  singly  supported  all  around  and  simply  supported-free  boumlary  condition. 

It  was  found  also  that  the  effect  of  through  the  thickness  shear  was  more  observable  at 
the  shallowo’  end  of  the  OMiiposite  shell  where  the  panel  bectxnes  thick  (S//i  <  1),  but 

u 

this  effect  disappears  as  the  curvature  effects  becomes  predominant 

Results  show  that  shallow  shell  thecny  is  a  good  approximation  when  the  S/c  ratio  is  less 
than  0.25,  but  as  the  ratio  increases  the  interaction  of  the  in-plane  displacements  must  be 
considered  as  part  of  the  eigenvector-eigenvalue  solution. 
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Appemlix  A:  More  Dctaib  Related  to  Um  Subepace  IteratkNi  Method 


The  purpose  of  this  section  is  to  present  stmie  additional  details  in  the  subqiace  method 
used  in  die  DSHELL  program. 


A.1  Subapace  Iteratkm  Method 

The  subspace  iteration  algorithm  was  developed  by  Bathe  [1]  and  is  explained  in  this 
section.  The  goal  of  the  subspace  iteration  method  is  to  sdive  for  the  ”p"  lowest 
eigenvalues  and  eigenvectors  which  satisfy 

K^  =  M4>A  (A.l) 

whereO  =  [(Oi}  (  ^}....(^p)  ],  A=diagmat]ixofei^nvaluesXi's  (i  =  1, 2, 
and  K  and  M  are  the  symmetric  stiffness  and  mass  matrices  respectively.  The  starting 
iterationvectOTS,  x, ,  span  the  subspace  (eigensubspace)£,  and  iteration  continues  until, 

widi  sufficient  accuracy,  £.is  spanned.  £.  is  the  p-dimensional  subspace  spanned  by  the 
"p"  lowest  eigenvectors.  The  starting  iteration  vectors  should  be  chosen  to  excite  those 
degrees  of  freedom  possessing  large  mass  and  mall  stiffness  [1:682].  One  method  for 
obtaining  x^ ,  which  has  bem  effectively  used,  begins  by  letting  the  first  column  of  [M]jc, 
be  the  majev  diagmial  of  [M].  This  establishes  excitation  of  all  degrees  of  freedom 
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possessing  mass.  The  other  columns  [M]x,  are  unit  vecttxs  with  a  "1”  at  locations 
with  the  smallest  KoAny  ratios.  Therefore,  all  locations  but  one  within  these  columns 
contains  zeros.  Hie  last  column  in  can  be  a  random  vector.  After  the  starting 
iteration  vector  is  determined,  iteration  begins  from  to  £^is  the  eventual  desired 
result ),  where  A  =  1, 2, 3, ...,  n  by  applying 

ATi  =  Mx*  (A.2) 

As  used  here,  "n"  is  the  number  of  iterations  required  to  obtained  convergence.  Equation 
( A.1 )  is  solved  for  x  by  premultiplication  of  both  sides  of  the  equation  by  [AT]'*  to 

yield 

(A.3) 

Next,  the  projections  of  the  AT  and  M  operators  on  the  subspace  are  determined  by 

[1] 


K,.,  = 

■*  Ll  M 

(A.4) 

Af*..  = 

XmMx 

(A.5) 

By  using  the  above  relations  for  AT^^,  and  the  eigensystem  of  the  f  ejected  operators 
can  be  solved  by  [1] 

^  ft+i  Q  M  —  M  Q  A  ( A.6  ) 
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Equation  ( A.5 )  represents  the  reduced  eigenvalue  problem  fOT  the  subspace  and  can 
be  solved  by  any  popular  technique  (DSHELL  uses  Householder's  Method  followed  by  a 
LR  algiHithm)  to  obtained  values  fen:  Qm  and  An  improved  iq)pn>ximation  to  the 
eigenvectors  is  found  from  [1] 

•*  *+I  ~  M  Q  M  (  ) 

This  procedure  is  iterated  until  the  desired  convergence  is  obtained.  This  x  is  now 

substituted  into  Equation  (  A.2 )  in  place  of  and  the  algorithm  iterated  until  sufficient 
convergence  has  been  obtained.  HgureA-1  presents  a  flow  chart  of  the  subspace 

a 

iteration  algorithm  incorporated  in  DSHELL.  In  subspace  iteration,  the  iteration  vectors 
are  ordered  such  that  the  vectors  converging  to  stmed  as  the  first, 

second, ...,  p*  columns  of  x  ,  respectively.  As  long  as  no  vector  contained  within  x, 
isorthogonal  to  any  of  the  required  eigenvectms,  then  as  k  tends  to  infinity,  A;^^, 
approaches  to  A  ( the  required  diagonal  matrix  of  "p"  eigenvalues  )  and  x  approaches 
[<|>]  ( the  required  column  matrix  cS  "p"  ergenvectOTS )  [1] 
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FigA-1  Subspace  Iteration  Algorithm 
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